LOW MACH NUMBER LIMIT FOR THE MULTI-DIMENSIONAL 
FULL MAGNETOHYDRODYNAMIC EQUATIONS 
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Abstract. The low Mach number limit for the multi-dimensional full magnc- 
tohydrodynamic equations, in which the effect of thermal conduction is taken 
into account, is rigorously justified in the framework of classical solutions with 
small density and temperature variations. Moreover, we show that for suffi- 
ciently small Mach number, the compressible magnctohydrodynamic equations 
admit a smooth solution on the time interval where the smooth solution of the 
incompressible magnctohydrodynamic equations exists. In addition, the low 
Mach number limit for the ideal magnctohydrodynamic equations with small 
entropy variation is also investigated. The convergence rates are obtained in 
both cases. 



1. Introduction 

The magnetohydrodynamic (MHD) equations govern the motion of compressible 
quasi-neutrally ionized fluids under the influence of electromagnetic fields. The full 
three-dimensional compressible MHD equations read as (see, e.g., [12,15,22,23]) 

dtp + div(pu) = 0, (1.1) 

d t {pu) +div(pu® u) + Vp= -!-(V x H) x H + divf, (1.2) 

47T 



d t U - V x (u x H) = -V x (zA7 x H), divH = 0, (1.3) 

1 

47T 



dtS + div (u(f + p)) = -!-div((u x H) x H) 



+ divf— H x (V x H) + u* + kW). (1.4) 

\4-7T / 

Here x £ fl, and f2 is assumed to be the whole R 3 or the torus T 3 . The unknowns p 
denotes the density, u = (ui,it2,«3) £ R 3 the velocity, H = (iJi, iJ 2 , -H3) £ K 3 the 
magnetic field, and 9 the temperature, respectively; VP is the viscous stress tensor 
given by 

* = 2/iP(u) + Adivu I 3 



Date: January 20, 2013. 

2000 Mathematics Subject Classification. 76W05, 35B40. 

Key words and phrases. Full MHD equations, smooth solution, low Mach number limit. 

1 



2 



SONG JIANG, QIANGCHANG JU, AND FUCAI LI 



with D(u) = (Vu + Vu T )/2, I3 being the 3x3 identity matrix, and Vu T the 
transpose of the matrix Vu; £ is the total energy given by £ = £' + |H| 2 / (87r) and 
£' = p (e + |u| 2 /2) with e being the internal energy, p|u| 2 /2 the kinetic energy, 
and |H| 2 /(87r) the magnetic energy. The viscosity coefficients A and p of the flow 
satisfy 2/i + 3A > and p, > 0; v > is the magnetic diffusion coefficient of the 
magnetic field, and k > is the heat conductivity. For simplicity, we assume that 
p, A, v and k are constants. The equations of state p = p(p, 9) and e = e(p, 9) relate 
the pressure p and the internal energy e to the density p and the temperature 9 of 
the flow. 

The MHD equations have attracted a lot of attention of physicists and math- 
ematicians because of its physical importance, complexity, rich phenomena, and 
mathematical challenges, see, for example, [2, 4-6, 8, 9, 11—13, 23, 31] and the refer- 
ences cited therein. One of the important topics on the equations (1.1)-(1.4) is to 
study its low Mach number limit. For the isentropic MHD equations, the low Mach 
number limit has been rigorously proved in [14, 16, 17, 20]. Nevertheless, it is more 
significant and difficult to study the limit for the non-iscntropic models from both 
physical and mathematical points of view. 

The main purpose of this paper is to present the rigorous justification of the 
low Mach number limit for the full MHD equations (1.1)- (1.4) in the framework of 
classical solutions. 

Now, we rewrite the energy equation (1.4) in the form of the internal energy. 
Multiplying (1.2) by u and (1.3) by H/(47r), and summing them together, we 
obtain 

|(>l 2 + ^|HP) + idiv( H uPu) + V ! ,.u 

= div* ■ u + — ( V x H) x H • u + — V x (uxH)-H 

47T 47T 

— V x (V x H) • H. (1.5) 

47T 

Using the identities 

div(H x(VxH)) = |VxH| 2 -Vx(VxH)-H (1.6) 

and 

div((u x H) x H) = (V x H) x H • u + V x (u x H) • H, (1.7) 
and subtracting (1.5) from (1.4), we obtain the internal energy equation 

dt(pe) + div(pue) + (divu)p = -^|V x H| 2 + * : Vu + kA9, (1.8) 

47T 

where "J : Vu denotes the scalar product of two matrices: 

3 2 

* : Vu = E f + S 1 ) + A l divu ! 2 = 2m|B(u £ )| 2 + A(trB(u e )) 2 . 
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In this paper, we shall focus our study on the ionized fluid obeying the perfect 
gas relations 

p = <Rp9, e = c v 9, (1.9) 

where the constants 9t, cy > are the gas constant and the heat capacity at constant 
volume, respectively. We point out here that our analysis below can be applied to 
more general equations of state for p and e by employing minor modifications in 
arguments. 

To study the low Mach number limit of the system (1.1)— (1.3) and (1.8), we use 
its appropriate dimensionlcss form as follows (see the Appendix for the details) 



d t p + div(pu) = 0, (1.10) 

p(d t u + u ■ Vu) + Yil^l = (V x H) x H + div*, (1.11) 

3,H-Vx(uxH) = -Vx(i/VxH), divH = 0, (1.12) 

p(d t 9 + u • V0) + (7 - 1)/J0divu = e 2 t/|V x H| 2 + e 2 * : Vu + nA6, (1.13) 



where e = M is the Mach number and the coefficients p, A, v and k are the scaled 
parameters. 7=1 + 91/cy is the ratio of specific heats. Note that we have used 
the same notations and assumed that the coefficients p, A, v and k are independent 
of e for simplicity. Also, we have ignored the Cowling number in the equations 
(1.10)— (1.13), since it does not create any mathematical difficulties in our analysis. 

We shall study the limit as e — > of the solutions to (1.10)-(1.13). We further 
restrict ourselves to the small density and temperature variations, i.e. 

p=l + eq, e = l + ecj). (1.14) 

We first give a formal analysis. Putting (1.14) and (1.9) into the system (1.10)- 
(1.13), and using the identities 

curl curl H = V divH - AH, 
V(|H| 2 ) = 2H-VH + 2HxcurlH, (1.15) 
curl (u x H) = u(divH) - H(divu) + H Vu u VH, (1.16) 

then we can rewrite (1.10)-(1.13) as 

d t q e + u e • V<f + -(1 + eg e )divu e = 0, (1.17) 

(1 + eq")(d t u e + u £ • Vu e ) + - [(1 + e<Z e )V0 e + (1 + e<£ e )V<f ] 

- H £ ■ VH e + iv(|H e | 2 ) = 2^div(B(u e )) + AV(trD(u £ )), (1.18) 
d t H c + u e • VH e + divu e H £ - H e • Vu e = i/AH c , divH e = 0, (1.19) 
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(l + eg 6 )($0 e +u e - V0 £ ) + 



7 -l 



(l + e g £ )(l + e0 e )divu e 



r 



kA0 £ + e(2 M |B(u e )| 2 + A(trD(u e )) 2 ) + z/e|V x H e | 2 . 



(1.20) 



Here we have added the superscript e on the unknowns to stress the dependence 
of the parameter e. Therefore, the formal limit as e —> of (1.17)— (1.19) is the 
following incompressible MHD equations (suppose that the limits u e — > w and 
H £ -> B exist.) 



In this paper we shall establish the above limit rigorously. Moreover, we shall 
show that for sufficiently small Mach number, the compressible flows admit a 
smooth solution on the time interval where the smooth solution of the incompress- 
ible MHD equations exists. In addition, we shall also study the low Mach number 
limit of the ideal compressible MHD equations (namely, fi = \ = v = K = 
in (1.1)— (1.4)) for which small pressure and entropy variations are assumed. The 
convergence rates are obtained in both cases. 

We should point out here that it still remains to be an open problem to prove 
rigorously the low Mach number limit of the ideal or full non-isentropic MHD equa- 
tions with large temperature variations in the framework of classical solutions, even 
in the whole space case, although the corresponding problems for the non-isentropic 
Euler and the full Navier-Stokes equations were solved in the whole space [1,25] or 
the bounded domain in [27]. The reason is that the presence of the magnetic field 
and its interaction with hydrodynamic motion in the MHD flow of large oscilla- 
tion cause serious difficulties. Wc can not apply directly the techniques developed 
in [1,25,27] for the Euler and Navier-Stokes equations to obtain the uniform esti- 
mates for the solutions to the ideal or full non-isentropic MHD equations. In the 
present paper, however, we shall employ an alternative approach, which is based 
on the energy estimates for symmetrizable quasilinear hyperbolic-parabolic systems 
and the convergence-stability lemma for singular limit problems [3,30], to deal with 
the ideal or full non-isentropic MHD equations. There are two advantages of this 
approach: The first one is that we can rigorously prove the incompressible limit 
in the time interval where the limiting system admits a smooth solution. The 
second one is that the estimates we obtained do not depend on the viscosity and 
thermodynamic coefficients, compared with the results in [10] where all-time exis- 
tence of smooth solutions to the full Navier-Stokes equations was discussed and the 
estimates depend on the parameters intimately. 



<9 f w + w ■ Vw + Vtt + -V(|B| 2 ) - B VB = /xAw, 
<9 f B + w • VB - B ■ Vw = i/AB, 



(1.21) 



(1.22) 
(1.23) 



divw = 0, divB = 0. 
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For large entropy variation and general initial data, the authors have rigorously 
proved the low Mach number limit of the non-isentropic MHD equations with zero 
magnetic diffusivity in [18] by adapting and modifying the approach developed 
in [25]. We mention that the coupled singular limit problem for the full MHD 
equations in the framework of the so-called variational solutions were studied re- 
cently in [21,26]. 

Before ending the introduction, we give the notations used throughout the cur- 
rent paper. We use the letter C to denote various positive constants independent of 
e. For convenience, we denote by H l = H l (il) (I <E R) the standard Sobolev spaces 
and write || • ||; for the standard norm of H l and || • || for || • 1 1 o - 

This paper is organized as follows. In Section 2 we state our main results. The 
proof for the full MHD equations and the ideal MHD equations is presented in 
Section 3 and Section 4, respectively. Finally, an appendix is given to derive briefly 
the dimensionless form of the full compressible MHD equations. 



We first recall the local existence of strong solutions to the incompressible MHD 
equations (1.21)— (1.23) in the domain J7. The proof can be found in [7,28]. Recall 
here that fi = R 3 or Q, = T 3 . 

Proposition 2.1 ( [7,28]). Let s > 3/2 + 2. Assume that the initial data (w, B)| t= o 
= (wo,Bo) satisfy Wq £ iJ s ,Bo € H s , and divwo = 0, divBo = 0. Then, there 
exist a T* € (0,oo] and a unique solution (w, B) G L°°(0, T*; H s ) to the incom- 
pressible MHD equations (1.21)— (1.23), and for any < T < T* , 



Denoting U £ = (q e , u c , H c , (j> e ) T , we rewrite the system (1.17)-(1.19) in the vec- 



2. Main results 



sup {||(w,B)(*)llff- + \\(d t w,d t B)(t)\\ Hs - 2 + ||V7r(t)|U.- 2 } < C. 



0<t<T 



tor form 



3 




(2.1) 



where 



Q(U e ) = (0, Ffu'j^Aff, K A(/) e + e(L(u e ) + G(H e ))) 



with 



F(u £ ) = 2 A idiv(D(u e )) + AV(trD(u e )) 
L(u e ) = 2 A i|B(u £ )| 2 + A(trP(u e )) 2 , 
G(H e ) = v\V x H £ | 2 , 



(i 
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and the matrices Aj(U e ) (0 < j < 3) are given by 









= diag(l, 1- 


he g £ ,l + eg e ,l 


+ e g £ ,l,l,l,l + e<f), 






A,{U*) = 
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It is easy to see that the matrices Aj(U e ) (0 < j < 3) can be symmetrized by 
choosing 

i (C/ £ ) = diag((l + e^)(l + eg 6 )" 1 , 1, 1, 1, 1, 1, 1, [(7 - 1)(1 + e^)}~ 1 ). 

Moreover, for U £ £ Gi CC G with G being the state space for the system (2.1), 
Aq(U € ) is a positive definite symmetric matrix for sufficiently small e. 
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Assume that the initial data U e (x 7 0) = Uq(x) = {qo(x), u Q (x), Hg(x), </> (x)) T £ 
H s and Uq(x) £ Go, Go CC G. The main theorem of the present paper is the 
following. 

Theorem 2.2. Let s > 3/2 + 2. Suppose that the initial data Uq(x) satisfy 
C/ o £ (x)-(0,w o ( 2; ),Bo(x),0) T =0(e). 

s 

Let (w,B,7r) be a smooth solution to (1.21)— (1.23) obtained in Proposition 2.1. If 
(w,tt) £ C{[0,T*},H s+2 ) n C 1 ([0,T*],H») with T* > finite, then there exists a 
constant eo > such that, for all e < eo, the system (2.1) with initial data Uq(x) 
has a unique smooth solution U e (x 1 t) £ C([0, T*\, H s ). Moreover, there exists a 
positive constant K > 0, independent of e, such that, for all e < eo, 



sup 

te[o,T* 



(^,w,B,|rr 



< Ke. (2.2) 



Remark 2.1. From Theorem 2.2, we know that for sufficiently small e and well- 
prepared initial data, the full MHD equations (1.1)-(1.4) admits a unique smooth 
solution on the same time interval where the smooth solution of the incompressible 
MHD equations exists. Moreover, the solution can be approximated as shown in 
(2.2). 

Remark 2.2. We remark that the constant K in (2.2) is also independent of the 
coefficients /i, v and k. This is quite different from the results by Hagstorm and 
Loranz in [10], where the estimates do depend on /i intimately. 

Our approach is still valid for the ideal compressible MHD equations. However, 
we will give a particular analysis for the ideal model with more general pressure 
by using the entropy form of the energy equation rather than the thermal energy 
equation in (1.8). 

The ideal compressible MHD equations can be written as 

dtp + div( i ou) = 0, (2.3) 
dt(pu) + div (pu ® u) + Vp = — (V x H) x H, (2.4) 

47T 

<9,H - V x (u x H) = 0, divH = 0, (2.5) 

3t£ + div (u(f ' +p)) = 7^-div((u x H) x H) . (2.6) 
With the help of the Gibbs relation 

and the identity (1.7), the energy balance equation (2.6) is replaced by 



0dS = de+pd 



dt(pS) + div(pSu) =0, (2.7) 
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where S denotes the entropy. We reconsider the equation of state as a function of 
S and p, i.e. p = R{S,p) for some positive smooth function R defined for all S and 
p > 0, and satisfying ^ > 0. Then, by utilizing (2.3), (1.15) and (1.16), together 
with the constraint divH = 0, the system (2.3)-(2.5) and (2.7) can be written in 
the dimcnsionlcss form as follows (after applying the arguments similar to those in 
the Appendix): 



R(S e ,p e )(dtu e + u £ • Vu £ ) + -4- - H £ • VH £ + -V(|H £ | 2 ) = 0, (2.9) 



A{S\p e )(d t p t + u £ • Vp £ ) + divu £ = 0, (2.8) 

Vp £ 1 

-|--H e -vir + -< 

d t W + u £ • VH £ + divu £ H £ - H £ • Vu £ = 0, divH 6 = 0, (2.10) 

d t S e + u £ • VS* £ = 0, (2.11) 

where A(S^p') = R(s \ pC) a %/\ 

To study the low Mach number limit of the above system, we introduce the 
transformation 

p e (x, t) = pePfW, S e (x, t) = S + e® e (x, t), (2.12) 
where p and S_ are positive constants, to obtain that 

a(S + e6 £ , eq c ){d t q e + u e • \7q c ) + -divu e = 0, (2.13) 

e 

r(S + e6 £ , eg £ )(9 t u £ + u £ • Vu £ ) + -X/q" - H £ • VH £ + iv(|H £ | 2 ) = 0, (2.14) 

9 t H £ + u £ • VH £ + divu £ H £ - H e • Vu £ = 0, divH £ = 0, (2.15) 

<9 t 6 £ + u £ • V6 £ = 0, (2.16) 
where 

pe eqe dR{S e ,s) 



a(S e ,eq e )=A(S e ,pe eq )pe eq 
r £ (5 £ , eg £ ) 



R{S^pe«i e ) ds 
R{S\pe tq ') 



s—pe eq 



pe eqe 



Making use of the fact that curlV = and letting e — > in (2.13) and (2.14), 
we formally deduce that divv = and 

curl (r(S, 0)(d t v + v • Vv) - (V x J) x J) = 0, 

where we have supposed that the limits u £ — > v and H £ J exist. Thus, we can 
expect that the limiting system of (2.13)-(2.16) takes the form 

r(S, 0)(d t v + v • Vv) - (V x J) x J + VII = 0, (2.17) 

<9 t J + v-VJ-J-Vv = 0, (2.18) 

divv = 0, divJ = (2.19) 



LOW MACH NUMBER LIMIT FOR FULL MAGNETOHYDRODYNAMIC EQUATIONS 9 



for some function IT. 

In order to state our result, we first recall the local existence of strong solutions 
to the ideal incompressible MHD equations (2.17)-(2.19) in the domain fl. The 
proof can be found in [7,28]. 

Proposition 2.3 ([7,28]). Let s > 3/2 + 1. Assume that the initial data (v, J)|t=o 
= (vo,Jo) satisfy vo G H s ,3q G H s , and divvo = 0, divJo = 0. Then, there 
exist a T* G (0,oo] and a unique smooth solution (v, J) G L°°(0, T*; H s ) to the 
incompressible MHD equations (1.21) -(1.23), and for any < T < T* , 

sup {\\(v,j)(t)\\ H s + ikc^cUX*)!!^-! + ||vn(t)|| ff .-i} < c. 

0<t<T 

In the vector form, we arrive at, for V e = (q e , u e , H e , O e ) T , that 

3 

A Q (eQ e ,eq e )d t V e +^ [u)A {e& ,eq e ) + e^C, + B^H 6 )}^ 6 = 0, (2.20) 
where 

A (e9 e , af) = diag(a(5 e , eq e ), r(5 e , eg e ), r(S\ eq% r(5 e , eq% 1, 1, 1, 1), 

and is symmetric constant matrix, and Bj(H e ) is a symmetric matrix of H c . 
Assume that the initial data for the equations (2.20) satisfy 

V %x) = (%(x),u (x),Hl(x),Q (x)) T &H S , and Vff(z) G Go, G CC G 

with G being state space for (2.20). Thus, our result on the ideal compressible 
MHD equations reads as 

Theorem 2.4. Let s > 3/2 + 1. Suppose that the initial data V e (x) satisfy 

\/ o e (x)-(0,vo(x),J o ( a; ),0) T =0(e). 

Let (v, J,II) be a smooth solution to (2.17)-(2.19) obtained in Proposition 2.3. If 
(v,IT) G C([0,T*],iP +1 ) n C^IO,^], H s ) with T* > /iraie, i/iera there exists a 
constant t\ > suc/i i/iai, /or all e< e\, the system (2.20) with initial data V e (x) 
has a unique solution V e (x,t) G C([0,T*], H s ). Moreover, there exists a positive 
constant K\ > such that, for all e < e±, 

sup \\V e {-,t) - (eIT,v, J,en) T || s < K x e. (2.21) 

te[o,T„] 

3. Proof of Theorem 2.2 

This section is devoted to proving Theorem 2.2. First, following the proof of the 
local existence theory for the initial value problem of symmetrizable hyperbolic- 
parabolic systems by Volpert and Hudjaev in [29], we obtain that there exists a 
time interval [0,T] with T > 0, so that the system (2.1) with initial data Uq(x) 
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has a unique classical solution U e (x,t) £ C([0,T},H s ) and U e (x,t) £ G2 with 



to prove the uniform boundedness of the solutions. See also [19] for some relative 
results. 

Now, define 



T t = sup{T > : U c (x,t) £ C([0, T],H S ), U e (x, t) £ G 2 ,V(x,t) £ ft x [0,T]}. 



Note that T £ depends on e and may tend to zero as e goes to 0. 

To show that lim c , T f > 0, we shall make use of the convergence-stability lemma 
which was established in [3, 30] for hyperbolic systems of balance laws. It is also 
implied in [30] that a convergence-stability lemma can be formulated as a part of 
(local) existence theories for any evolution equations. For the hyperbolic-parabolic 
system (2.1), we have the following convergence-stability lemma. 

Lemma 3.1. Let s > 3/2 + 2. Suppose that Uq(x) £ G ,G CC G, and Uq(x) £ 
H s , and the following convergence assumption (A) holds. 

(A) There exists T* > and U e <G L°°(0, T*; H s ) for each e, satisfying 



such that for t £ [0, minjT*, T e }), 

sup\U e (x,t) - U e (x,t)\ = o(l), sup\\U e {x,t) - U e (x,t)\\ s = O(l), as e -> 0. 

x,t t 

Then, there exist an e > such that, for all e £ (0, e], it holds that 



To apply Lemma 3.1, we construct the approximation U e = (q e , v e , B c , (f> e ) T with 
q c = e7r/2, v e = w, B £ = B, and £ = eir/2, where (w, B, tt) is the classical solution 
to the system (1.21)-(1.23) obtained in Proposition 2.1. It is easy to verify that U € 
satisfies 



G2 CC G. We remark that the crucial step in the proof of local existence result is 



\J{U e (x,t)}ccG, 



T e > r*. 




- B e • VB e + -V(|B e | 2 ) = /iAv e + — Tr(w t + w • Vw + Vtt), (3.2) 
dt~B e + v e • VB e + divv e B e - B e • Vv f = ^AB e , divB e = 0, (3.3) 



(l + e 9e )(<9 t (/> e +v e -V0 e ) + 



7 



1 



(l + e (7£ )(l + e0 e )divv e 



c 




(3.4) 
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We rewrite the system (3.1)-(3.4) in the following vector form 



A (U e )dtU t + ^A j (U € )d j U< = S(U e ) + R, 
3=1 



(3.5) 



with S(U e ) = (0,/iAv e ,i/AB e ,0) T and 



R 



f (wt + w • Vtt) 
^-7t(wj + w • Vw + Vn) 



,2 + 4 



V 



71") (Wt + W • V7r) 




Due to the regularity assumptions on (w,7r) in Theorem 2.2, we have 

max \\R(t)\\ s < Ce. 

t£[0,T«] 

To prove Theorem 2.2, it suffices to prove the error estimate in (2.2) for t E 
[0, min{T*, T e }) thanks to Lemma 3.1. To this end, introducing 

E = U e -U £ and Aj{U) = 1 (U)A j (U), 

and using (2.1) and (3.5), we see that 

3 

Et + YtMV^E*, =(Aj(U e ) - Aj(U e ))U eXj + ^o 1 ([/ e )Q([/ e ) 
i=i 

-Ao\U e )(S(U e )+R). (3.6) 

For any multi-index a satisfying |a| < s, we take the operator D a to (3.6) to 
obtain 

3 

d t D a E + ^2Aj(U f )d Xi D a E = Pf + P 2 Q + Q a + R a (3.7) 



with 



p r = » l .,<'";^ /r/ - - ^"(^r)^^}. 

P« =^£)«{(^.(C/ e ) -AjWU^}, 

3=1 

Q a = D a {A^\W)Q(W) - AoHUe)S(U e )}, 
R a = D a {A^ 1 (U e )R}. 



Define 



l + e 9 e ' l + e<f' l + e<f' (7- 1)(1 + e^) . 
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and the canonical energy by 

||P|| 2 := / (A (U e )E,E)dx. 



Note that Ao(U e ) is a positive definite symmetric matrix and Ao(U e )Aj(U e ) is 
symmetric. Now, if we multiply (3.7) with Ao(U e ) and take the inner product 
between the resulting system and D a E, we arrive at 

^-\\D a E\\l =2 / (TD a E,D a E)dx 
at J 

+ 2 J (D a E) T A (U e )(P? + P? + Q a + R a ), (3.8) 

where 

r = (d t , v) • (A ,A (u e )Ai(u*),A (u*)A 2 (u e ),A (u*)A 3 (u e )). 

Next, we estimate various terms on the right-hand side of (3.8). Note that our 
estimates only need to be done for t £ [0, min{T*, T e }), in which both U e and U € 
are regular enough and take values in a convex compact subset of the state space. 
Thus, we have 

C" 1 J \D a E\ 2 < \\D a E\\ 2 c <C J \D a E\ 2 (3.9) 

and 

\{D a E) T A^U^P? + P 2 Q + R a )\ < C(\D a E\ 2 + \P?\ 2 + |P 2 Q | 2 + |iT | 2 ). 

To estimate T, we write Aj(U e ) = u^I s + Aj{U t ). Notice that Aj{U e ) depends 
only on q e ,(j> e and H £ . Thus using (1.17) and (1.19), we have 

\T\ = — A + u) • Vi + i divu e + div(ioA(f/ £ )) 

= \Ao divu £ - A' 0m (1 + eg £ )divu £ - A' 0r)2 [(1 + e0 £ )divu £ + k(1 + eq 6 )' 1 A<f 

+ e 2 (L(u £ ) + G(H £ ))] + div(i ^-(L/ £ ))| 
<C + C7(|Vu £ | + |V<f| + |V0 £ | + |VH £ | + |A0 £ | + |Vu £ | 2 + |VH £ | 2 ) 
<C + C(\VE\ + \VE\ 2 ) + C\A((^ e -<f> c )\ + C{\VU £ \ + |Vf/ £ | 2 ) 

<c + c(||p|| s + ||p!| 2 ), 

where we have used Sobolev's embedding theorem and the fact that s > 3/2 + 2, 
and the symbols A 0m and A denote the differentiation of Aq with respect to p e 
and 9 e , respectively. 
Since 

A j {U*)d Xi D a E-D a {A j (U')d Xj E) = - ]T (fjd^AjiU^d^E^ 



0</3<a 
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0</3<a V 3 ' 

we obtain, with the help of the Moser-type calculus inequalities in Sobolev spaces, 
that 

||Pf|| <C{(l + ||(<H e )||.)||B a J| w _ 1 + ||e-H^/(ff e ,0 e )a a -^)||} 

+ C\\dP[(l + eq*)- 1 ^ - B £ ) + ((1 + eg 6 )" 1 - (1 + eq^B^d^ E Xj \\ 

<c(i + 1|^||. + ik^.^iidii^Hh-! 

<C(l + \\E\\ s s )\\E Xj \\ lal , 

where f{q e , </> £ ) = (1 + eg 6 ) + (7 - 1)(1 + e<t> c ) + (1 + etf)~ l (l + e<f). 

Similarly, utilizing the boundedness of ||C/ e || s+ i, the term P£ can be bounded as 
follows: 

\\m<C\\U tXl \\.\\A i {V*)-A i (V<)\\\ a \ 

<C\\(u^~v e3 )l 8 +A 3 (U*)-A 3 (U e )\\ lal 

< C(l + ||u e - v e || w + ||H 6 - B e || w ) + C\\e-\f {q'rf) - /(g e ,0e))|||„| 

< C(l + \\q e + V3 (q' - q e ) + <f> £ + r U (^ - £ )O|£|| W 

^ca+ii^ii^iiEii^i, 

where < 773,774 < 1 are constants. 

The estimate of J (D a E) T A (U e )Q a is more complex and delicate. First, we 
can rewrite J(D a E) T A (U e )Q a as 

J (D a E) T A (U e )Q a = J D a (u € - v e )D a [(1 + eq^F^) - fi(l + eq^Av.] 
+ v j D a (H. c — B c )(l + eq e )~ 1 D a (AH e — AB e ) 

+ «(7 - I)" 1 / - 4>e)(l + e^)- 1 D a {{l + eg)"^ - (1 + eqA^Afc} 

+ c(7 - I)" 1 / £> Q (0 e - + e0 e )- 1 D a {(l + eg 6 )- 1 ^^ 6 ) + G(H e ))} 



By integration by parts, the Cauchy and Moser-type inequalities, and Sobolev's 
embedding theorem, we find that Q u can be controlled as follows: 

Qu = J D a (u € - v e )D Q {(l + egTVA(u e _ Ve ) + ((U + A )Vdiv(u e - v e )} 
+ nj D a (tf - v e )D a {[(l + efy 1 - (1 + eg^-^AvJ 

<-J ^I^V^-Ve)! 2 -/ ^|^div ( U £ -V £ )P 
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+ /V(u e -v e ) D< > [(l + eq e )- 1 ]D a - f '{nA(u e -v e ) 

Q<P<a 

+ (/i + A)Vdiv(u e - v £ )} + C||£7||f„| + C\\E\\j + Ce\\D a W(u e - v e )|| 2 
<-C J M|£"V(u e -v e )| 2 -C yO* + A)|D a div(u e -v e )| a + q|S||f a | 

+ C7e||D a V(u e - v e )|| 2 + C\\E\\t + CWEWlWEW^ + j D a (u< v e )- 

D P [Q- + etf)- 1 ]B a - /) {/xA(u e - v e ) + (/x + A)Vdiv(u e - v £ )} 

Kj8<a 

< - C J ^|zrv(u e - v e )| 2 - C f(fx+ A)|D Q div(u e - v e )| 2 
+ Ce\\D a V(u e - v e )|| 2 + C\\E\\t + C\\E\\f al . 
Similarly, the terms Qh, Qfa and Q^ 2 can be bounded as follows: 
Qh < -Cv J \D a V{W - B e )| 2 + Ce\\D a V{W - B e )|| 2 + C\\E\\* + C\\E 
Q<t>i < -Ck J \D a V{^ - £ )| 2 + Ce\\D a V{tf - ^ e )|| 2 + C\\E\\t + C\\E\\fa 



2 

\a\> 



and 



Qfc < Ce\\D a V{^ - e )|| 2 + C\\E\\i + C\\E\tf a{ . 



Putting all the above estimates into (3.8) and taking e small enough, we obtain 
that 

j t \\D a E\\ 2 e +£ J |£ Q V(u e - v e )| 2 + ^ J |£ Q V(H e -B e )| 2 

+ k J \D a V(^ - </> e )| 2 < C\\R a \\ 2 + C(l + \\E\\l s )\\E\\U + \\Ef s , (3.10) 

where we have used the following estimate 

H j |D Q V(u e - v e )| 2 + (a* + A) j |D Q div(u e - v e )| 2 > £ J \D a V{u e - v e )| 2 

for some positive constant £ > 0. 

Using (3.9), we integrate the inequality (3.10) over (0,i) with t < min{T e ,T*} 
to obtain 



\\D a E(t)\\ 2 <C\\D a E(0)\\ 2 + C / \\R a (T)\\ 2 dT 

Jo 

+ C f {(l + \\E\\f)\\E\\l l + \\E\\i}(r)dr. 



Summing up this inequality for all a with \a\ < s, we get 

\\E(t)\\ 2 s <C\\E(0)f s +C j T \\R(r)\\ldr + C f {(1 + \\E\\f)\\E\\l}(r)dr. 
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With the help of Gronwall's lemma and the fact that 

11^(0)112+ f \\R(t)f s dt = 0(e 2 ), 
Jo 

we conclude that 

\\E(t)f s < Ce'cxp |(7 + ||25(r)||*)dr} = S(f). 

It is easy to see that <J>(i) satisfies 

$'(t) = C(l + \\E(t)\\ 2 s s )$(t) < C$(t) + C$ s+1 (i). 

Thus, employing the nonlinear Gronwall-type inequality, we conclude that there 
exists a constant K, independent of e, such that 

\\E(t)\\ s < Ke, 

for all t E [0,min{T e ,T*}), provided $(0) = Ce 2 < exp(-CT*). Thus, the proof is 
completed. 

4. Proof of Theorem 2.4 

The proof of Theorem 2.4 is essentially similar to that of Theorem 2.2, and we 
only give some explanations here. The local existence of classical solution to the 
system (2.20) is given by the proof of Theorem 2.1 in [24]. For each fixed e, we 
assume that the maximal time interval of existence is [0,T e ). To prove Theorem 
2.4, it is crucial to obtain the error estimates in (2.21). For this purpose, we 
construct the approximation V e = (q e , v c , J e , 9 e ) T with q e = ell,v e = v,J e = J, 
and e = dl. It is then easy to verify that V e satisfies 

a(S_ + e<d e ,eq e )(d t qe + v e • Vg £ ) + -divv, 

e 

= ea(5 + e 2 n,e 2 II)(IT +v VII), (4.1) 
r(S + e9 e , eq t )(d t w t + v e ■ Vv £ ) + -Wq, - J e • VJ e + i(|J £ | 2 ) 

= [r(5 + ee e> e5 £ )-r(S,0)](v t + v.Vv) > (4.2) 
d t J e + v e • VJ e + divv e J £ - J e • Vv e = 0, divJ £ = 0, (4.3) 

9te e + v e -ve e = e(n t + v-vn). (4.4) 

Thus we can rewrite (4.1)-(4.4) in the vector form of (2.20) with a source term. 
Letting E = V e — V e , we can perform the energy estimates similar to those in the 
proof of Theorem 2.2 to show Theorem 2.4. Here we omit the details of the proof 
for conciseness. 
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5. Appendix 

We give a dimensionless form of the system (1.1)-(1.3) and (1.8) for the ionized 
fluid obeying the perfect gas relations (1.9) by following the spirit of [12]. Introduce 
the new dimensionless quantities: 

x t u 

H, = Jt *-A«.-f 

no po Oq 

where the subscripts denote the corresponding typical values and * denotes dimen- 
sionless quantities. For convenience, all the coefficients are assumed to be constants. 
Thus, the dimensionless form of the system (1.1)— (1.3) and (1.8) is obtained by a 
direct computation: 

^ + div*(/J*u*) = 0, 

P*^ + ]j^V*(M*) = C(V* x H*) x H* + ^div***, 

P*^T + (7 - lKMiv*u* = l2_l) C i\/ 2 |V, x H*| 2 
+ (7~1)M 2 , 

R RP r ' 

OH. 1 

— -i - V* x (u* x H*) = —V* x (V* x H*), div*H* = 0, 

where we have used the material derivative 

d d 

-jr = tt— + u* • V*, 

at* ot* 

and the new viscous stress tensor 

= 2B*(u*) + -div^u* I 3 

with B*(u*) = (V*u* + V*uJ)/2. 

In the above dimensionless system, there are following dimensionless character- 
istic parameters: 

Reynolds number: R = ^° — -, Mach number: M = — , 

fx a 

Prandtl number: P r = 2eEl^ magnetic Reynolds number: R m = — — -, 
k v 

Cowling number: C — ^ , 

where c p is the specific heat at constant pressure and do = \/V\9q is the sound 
speed. Note that 91 = c p — cy and 7 = c p /cy. 
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